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In this talk we discuss a purely numerical approach to next-to-leading order calculations in QCD. We present 
a simple formula, which provides a local infrared subtraction term for the integrand of a one-loop amplitude. In 
addition we briefly comment on local ultraviolet subtraction terms and on the required deformation of the contour 
of integration. 



1. Introduction 

The experiments at the LHC are faced with a 
high QCD jet rate. An accurate description of jet 
physics is therefore mandatory. Although jet ob- 
servables can rather easily be modelled at leading 
order (LO) in perturbation theory, this descrip- 
tion suffers several drawbacks. A leading order 
calculation depends strongly on the renormalisa- 
tion scale and can therefore give only an order-of- 
magnitude estimate on absolute rates. Secondly, 
at leading order a jet is modelled by a single par- 
ton. This is a very crude approximation and over- 
simplifies inter- and intra-jet correlations. The 
situation is improved by including higher-order 
corrections in perturbation theory. Of particu- 
lar phenomenological interest are next-to-leading 
order (NLO) corrections to multi-jet observables, 
where the number of jets is in the range of 3 to 6 
or 7. This is not a simple task: The complexity of 
the calculation increases with the number of final 
state jets. There has been significant progress 
in the last years with these calculations, either 
based on traditional Feynman diagram techniques 
[1I2I3I4I5I6I7I8| or on cut techniques |9I10I11|I^ . 

In this talk we discuss a third method which 
is based on a numerical integration for the loop 
integrals. This by itself is not a new idea and has 
been discussed before [13 14 15 16 17 18,. What 
is new is a simple and compact formula which 
approximates the integrand of a one-loop QCD 
amplitude in all infrared singular regions. This 



approximation can be taken as a local countert- 
erm and opens the door for an efficient numerical 
implementation. In addition we need local coun- 
terterms for the ultraviolet divergences and a de- 
formation of the integration contour. We briefly 
comment on these points. 

2. The subtraction method 



The subtraction method |19I20I21] is widely 
used to render the real emission part of an NLO 
calculation suitable for a numerical Monte Carlo 
integration. At NLO one has the following con- 
tributions: 
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Here a rather condensed notation is used, da^ de- 
notes the real emission contribution, whose ma- 
trix element is given by the square of the Born 
amplitudes with {n + 3) partons. da^ gives 
the virtual contribution, whose matrix element is 
given by the interference term of the one-loop am- 
plitude with {n + 2) partons with the correspond- 
ing Born amplitude, da*^ denotes a coUinear sub- 
traction term, which subtracts the initial-state 
coUinear singularities. Taken separately, the indi- 
vidual contributions are divergent and only their 
sum is finite. In order to render the individual 
contributions finite, such that the phase space 
integrations can be performed by Monte Carlo 
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methods, one adds and subtracts a suitable cho- 
sen piece: 
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The term {On+idcr^ — Ondcr^) in the first bracket 
is by construction integrable over the (n + 1)- 
particle phase space and can be evaluated numer- 
ically. The subtraction term can be integrated 
over the unresolved one-parton phase space. Due 
to this integration all spin-correlations average 
out, but colour correlations still remain. In a 
compact notation, the result of this integration 
is often written as 



da"- 



da' 



I ® da^ -f K (g) da" + 'P ®da 



B 



(3) 



The notation (g) indicates that colour correlations 
due to the colour charge operators still re- 
main. The terms with the insertion operators K 
and P pose no problem for a numerical evalua- 
tion. The term I ® da^ lives on the phase space 
of the n-parton configuration and has the appro- 
priate singularity structure to cancel the infrared 
divergences coming from the one-loop amplitude. 
Therefore da^ da^ is infrared finite. 

The virtual correction term da^ contains the 
one-loop amplitude and is usually considered to 
be the bottleneck of a NLO calculation, da^ is 
given by 



da^ = 2 Re 
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A^^'^ denotes the renormalised one-loop ampli- 
tude. It is related to the bare amplitude by 
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,4p^ denotes the ultraviolet counterterm from 
renormalisation. The bare one-loop amplitude in- 
volves the loop integration 
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where denotes the integrand of the bare one- 
loop amplitude. Within our approach we extend 
the subtraction method to the integration over 
the virtual loop momentum k. We rewrite eq. (O 
as 
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The subtraction terms A^^Jf^, -^con ^"^d A^y are 
chosen such that they match locally the singu- 
lar behaviour of the integrand of ^[^l-g in D di- 
mensions. The expression in the first bracket can 
therefore be integrated numerically over the loop 
momentum k in four dimensions. The term A[}]ft 

so J z 

approximates the soft singularities, A'^^Jh approx- 
imates the coUinear singularities and the term 
A^y approximates the ultraviolet singularities. 
These subtraction terms have a local form simi- 
lar to eq. ©: 

d^k 

' gi^\ X ^ soft, coll, uv. (8) 
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The building blocks of the subtraction terms are 
process- independent. When adding them back, 
we integrate analytically over the loop momen- 
tum k. The result can be written as 



2 RjG ^ ^"^OT' 
= L(E)da^. 
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The insertion operator L contains the explicit 
poles in the dimensional regularisation parameter 
related to the infrared singularities of the one- 
loop amplitude. These poles cancel when com- 
bined with the insertion operator I: 



(I 4- L) «) da^ = finite. 
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The operator L contains, as does the operator I, 
colour correlations due to soft gluons. 

3. Colour decomposition 

It is convenient to decompose a full one-loop 
QCD amplitude into primitive amplitudes: 
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The colour structures are denoted by Cj while 

the primitive amplitudes are denoted by ^j-^'- In 
the colour-flow basis |22I23I24] the colour struc- 
tures are linear combinations of monomials in 
Kronecker (Sy 's. Primitive amplitudes are defined 
as a colour-stripped gauge-invariant set of Feyn- 
man diagrams with a fixed cyclic ordering of the 
external partons and a definite routing of the ex- 
ternal fermion lines through the diagram |25| . 

It is simpler to work with primitive one-loop 
amplitudes instead of a full one-loop amplitude. 
Our method exploits the fact that primitive one- 
loop amplitudes have a fixed cyclic ordering of the 
external legs and that they are gauge-invariant. 
The first point ensures that there are at maxi- 
mum n different loop propagators in the problem, 
where n is the number of external legs, while the 
second property of gauge invariance is crucial for 
the proof of the method. We therefore consider in 
the following just a single primitive one-loop am- 
plitude, which we denote by , while keeping in 
mind that the full one-loop amplitude is just the 
sum of several primitive amplitudes multiplied by 
colour structures. 

We introduce some notation related to the 
primitive amplitude A^^\ Since the cyclic order- 
ing of the external partons is fixed, there are only 
n different propagators occurring in the loop in- 
tegral. We label the external momenta clockwise 
by pi, p2, Pn and define q., = pi +p2 + ■■■+Pi, 
ki = k — Qi. We can write the bare primitive 
one-loop amplitude in Feynman gauge as 
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^bare integrand of the bare one-loop am- 

plitude. P{k) is a polynomial in the loop mo- 
mentum k. The prescription instructs us to 
deform - if possible - the integration contour into 
the complex plane to avoid the poles at kf = m^. 
If a deformation close to a pole is not possible, 
we say that the contour is pinched. If we restrict 
ourselves to non-exceptional external momenta, 
then the divergences of the one-loop amplitude 
related to a pinched contour are either due to 



soft or collinear partons in the loop. These diver- 
gences are regulated within dimensional regulari- 
sation by setting the number of space-time dimen- 
sions equal to D = 4 — 2e. A primitive amplitude 
which has soft or collinear divergences must have 
at least one loop propagator which corresponds 
to a gluon. An amplitude which just consists of 
a closed fermion loop does not have any infrared 
divergences. We denote by Ig the set of indices i, 
for which the propagator i in the loop corresponds 
to a gluon. 

4. Local infrared subtraction terms 

We can now present a formula for the local sub- 
traction terms related to the infrared divergences 
|26j . For simplicity we present here the formula 
for massless QCD. The extension to massive par- 
ticles is also known. The infrared subtraction 
term can be written in unintegrated form as a 
soft and a collinear part: 
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The Born partial amplitude Af''^ depends on the 
external momenta, but not on the loop momen- 
tum. The constants Si are given by 5^ = S*^ = 1 
and Sg ~ 1/2, the index i of Si refers to the exter- 
nal particles. The function guv ensures a regular 
behaviour of the collinear term in the ultraviolet 
region. A possible choice is [TBI 



guv {kf_i,k'^,fil) = i 



/ic is an arbitrary scale. Integrating the soft and 
the collinear part we obtain 
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Se = (47r)'^e~^'''^ is the typical volume factor of 
dimensional regularisation, 7^ is Euler's constant 
and /X is the renormalisation scale. The conditions 
under which a single diagram leads to an infrared 
divergence are well known [27128129] . The proof 
of eq. given in f26] uses then the fixed cyclic 
ordering and gauge invariance to establish the re- 
lation at the level of primitive amplitudes. 

Equation (|13p approximates the integrand of 
a primitive one-loop QCD amplitude in all soft 
and coUinear limits. The approximation is given 
by simple scalar two- and three-point functions, 
multiplied by a Born partial amplitude. One eas- 
ily observes that the integrated form in eq. (fT4|) 
agrees in the pole terms with the known result for 
the pole terms of a primitive one-loop amplitude. 

Equation (jl3p is a significant improvement with 
respect to the subtraction terms given in ref. |15) . 
In ref. |15j the subtraction terms where given 
graph by graph. In contrast to this we see that 
eq. ([T5|) is formulated at the level of amplitudes 
and not at the level of graphs. The improvement 
results from the fact that an amplitude can be cal- 
culated efficiently using recurrence relations. The 
recursive method is significantly faster compared 
to an approach based on individual diagrams. 

A few remarks are in order: Contrary to the 
subtraction terms for the real emission, there are 
no spin correlations in the subtraction terms for 
the integrand of the one-loop amplitude. In the 
real emission case spin correlations occur in the 
coUinear limit. This can be understood as fol- 
lows: From the proof of eq. one sees that 
in the coUinear limit always one of the coUinear 
gluons carries an unphysical longitudinal polari- 
sation. Hence, there are no correlations between 
two transverse polarisations of a gluon. 

Furthermore and also contrary to the subtrac- 



tion terms for the real emission, there is no depen- 
dence on the variant of dimensional regularisa- 
tion (conventional dimensional regularisation, 't 
Hooft-Veltman scheme, four-dimensional scheme) 
in the integrated result. At first sight this is puz- 
zling: In the real emission case the variant of 
dimensional regularisation introduces a scheme- 
dependent finite term. On the other hand uni- 
tarity requires that this scheme-dependent finite 
term cancels in the final result. It is instruc- 
tive to investigate how this cancellation occurs. 
The solution comes from the LSZ-reduction for- 
mula: The renormalised one-loop amplitude with 
fig gluons, Uq quarks and nq antiquarks is related 
to the bare amplitude by 

X Aare (pi , ■■■,Pn, Z^S^V^^^s) • (15) 

Zg is the renormalisation constant for the strong 
coupling, given by 

z, = , + 2.(_i|)l,c,„;,. (le) 

Z2 is the quark field renormalisation constant and 
Z3 is the gluon field renormalisation constant. 
The LSZ reduction formula instructs us to take 
as field renormalisation constants the residue of 
the propagators at the pole. In dimensional reg- 
ularisation this residue is 1 for massless parti- 
cles and therefore the field renormalisation con- 
stants are often omitted from eq. (|15p. However 
■2^2 = -^3 = 1 is due to a cancellation between ul- 
traviolet and infrared divergences |30] . In Feyn- 
man gauge we have 

Z2 = l + ^cJ- l-)+0{al), 

^3 = 1 + ^{2Ca-m{- 

47r \£iR Euv) 

+0(a2). (17) 

Here we indicated explicitly the origin of the 1/e- 
poles. These poles introduce scheme-dependent 
finite terms of ultraviolet and infrared origin 
with opposite sign. Now the cancellation of the 
scheme-dependent finite term of infrared origin 
is as follows: The scheme-dependent finite term 
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from the real emission contribution cancels with 
the scheme-dependent finite term of infrared ori- 
gin from the renormalisation constants. This 
leaves the scheme-dependent finite terms of ul- 
traviolet origin in the renormalisation constants 
and in the bare one-loop amplitude. These re- 
main and give the result of the calculation in the 
chosen scheme. One can convert from one scheme 
to another by a finite renormalisation. We remark 
that the scheme-dependence of the bare one-loop 
amplitude is entirely of ultraviolet origin. 

5. Local ultraviolet subtraction terms 

A primitive one-loop amplitude contains, apart 
from infrared divergences, also ultraviolet diver- 
gences. In an analytic calculation regulated by di- 
mensional regularisation these divergences man- 
ifest themselves in a single pole in the dimen- 
sional regularisation parameter e. If we would 
do the calculation within a cut-off regularisation 
we would find a logarithmic dependence on the 
cut-off. 

Within the numerical approach proposed here 
we need subtraction terms, which approximate 
the ultraviolet behaviour of the integrand locally. 
We first note that the one-loop amplitude has in 
a fixed direction in loop momentum space up to 
quadratic UV-divergences. These are reduced to 
a logarithmic UV-divergence only after angular 
integration. For a local subtraction term we have 
to match the quadratic, linear and logarithmic di- 
vergence. The ultraviolet subtraction terms have 
the form of propagator and vertex counterterms. 
(Diagrams with five or more loop propagators are 
never UV-divergent.) We can choose them such 
that after integration they are proportional to 
their tree- level counterpart. This ensures that 
the sum of all integrated UV subtraction terms 
is again proportional to a tree amplitude. As 
an example we give here the ultraviolet subtrac- 
tion term for the quark-photon vertex. In unin- 
tegrated form this term reads [31^ 

Vqq^ = iQ.eg^S^^fi^-'' (18) 

d^k 2 (1 - e) ^7^p -I- 4^2^v'7'' 



fiuv is an arbitrary scale. The term proportional 



to 4:^fjyj'^ in the numerator is not divergent, but 
ensures that the integrated expression has a sim- 
ple form. Integration leads to 
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+ Oie). (19) 



Similar expressions can be found for all other UV- 
divergent self-energy and vertex corrections. A 
possible choice of UV-subtraction terms has been 
given in [15] . 

6. Contour deformation 

Having a complete list of ultraviolet and in- 
frared subtraction terms at hand, we can en- 
sure that the integration over the loop momen- 
tum gives a finite result and can therefore be per- 
formed in four dimensions. However, this does 
not yet imply that we can simply or safely in- 
tegrate each of the four components of the loop 
momentum fc'^ from minus infinity to plus infinity 
along the real axis. There is still the possibility 
that some of the loop propagators go on-shell for 
real values of the loop momentum. If the contour 
is not pinched this is harmless, as we may escape 
into the complex plane in a direction indicated by 
Feynman's -l-iJ-prescription. However, it implies 
that the integration should be done over a region 
of real dimension 4 in the complex space C^. 

Let us consider an integral corresponding to a 
primitive one-loop amplitude with n propagators 
minus the appropriate IR- and UV-subtraction 
terms: 
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P{k') is a polynomial or more general a well- 
behaved function in the loop momentum k' . The 
integration is over a complex contour in order to 
avoid " whenever possible - the poles of the prop- 
agators. We set 

k'^ = k^' + iK''{k), (21) 
where fc'^ is real. This introduces a Jacobian 
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and the integral becomes 

J (27r) {k[{k)) -mf 

The efficiency of the Monte Carlo integration de- 
pends crucially on the function which de- 
fines the deformation. Possible choices have been 
discussed in [16|17] . 

7. Conclusions 

In this talk we discussed a numerical approach 
to NLO calculations in QCD. Our main new re- 
sult is a simple and compact formula, which ap- 
proximates the integrand of a primitive one-loop 
QCD amplitude in all limits where the loop inte- 
gration leads to a soft or collinear singularity. It 
thus serves as a subtraction term for the infrared 
divergences of a one-loop amplitude. When com- 
bined with the corresponding subtraction terms 
for the ultraviolet divergences we may perform 
the loop integration of the subtracted integrand 
in four dimensions. It should be mentioned that 
there is still the need to deform the integration 
contour into the complex plane, since we have 
to avoid singularities of the integrand in regions 
where the integration contour is not pinched. 
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